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THERMAL RADIATIVE PROPERTY MEASUREMENTS USING 
CYCLIC RADIATION 
by John R. Jack and Ernie W. Spisz 
Lewis Research Center 

SUMMARY 

An experimental technique that determines simultaneously the total hemispherical 
emittance, solar absorptance, and the absorptance-emittance ratio for various materials 
over a substantial temperature range is discussed. The method is based upon comparing 
the temperature response of a thin, flat, isothermal sample to a sinusoidal perturbation 
of an imposed incident radiation. An experimental validation of the method is presented 
along with experimental thermal radiative property data for a variety of materials at 
temperatures ranging from 150 to 600 K. Twelve metals, a Z-93 coating, a black velvet 
paint, and a composite material (a CdS solar cell) were investigated to demonstrate the 
capabilities and the general utility and versatility of the cyclic radiation method for de- 
termining thermal radiative properties. 


INTRODUCTION 

The design of space vehicles and space equipment requires accurate total radiative 
property data on all types of materials over a wide range of environmental conditions. 
This need stems from the fact that passive spacecraft temperature control is generally 
achieved by a radiation heat exchange between spacecraft components and the space en- 
vironment. Neither existing experimental data nor theoretical predictions are adequate 
to fill this need. Therefore, the designers responsible for spacecraft thermal control 
must rely completely upon accurate experimental determinations of thermal radiative 
properties that satisfy the required combination of thermal and environmental boundary 
conditions. 

Various techniques are available for the determination of thermal radiative proper- 
ties; most of which require a separate experiment for each radiative property desired. 
Of the many methods available, all can be classified as either a thermal or an optical 


technique. However, each method has associated with it experimental difficulties and 
limitations. 

Thermal methods are generally calorimetric and determine total radiative proper- 
ties either by equilibrium or transient temperature measurements. The calorimetric 
approach has inherent errors, such as accurate temperature measurement, heat leaks, 
background radiation, and convective heat-transfer effects; all of which can significantly 
affect the experimental results (ref. 1). This is especially true at low temperature 
levels, where the heat transferred by radiation is small. 

With optical methods, the spectral reflectance is measured over the appropriate 
spectral range and the desired thermal radiative property calculated. Such calculations 
are lengthy and arduous, but more importantly, the spectral reflectance method involves 
a series of assumptions related to the equipment and sample characteristics such as 
specularity, angular dependency, and polarization (ref. 2). In addition, although the re- 
flectance can be measured accurately, the error in the calculated thermal radiative 
properties a, e can be very large for highly reflecting materials. For example, pol- 
ished aluminum with a reflectance of 0. 98 ± 0. 005 at a wavelength of 10 micrometers 
yields an error in spectral emittance of 25 percent. 

Thus, although some information on the thermal radiation characteristics of mate- 
rials is available, the data are, in general, not completely reproducible from specimen 
to specimen nor from source to source for the same material. This, coupled with the 
facts that our basic understanding of the physical phenomena is not sufficient to permit 
reliable theoretical predictions and that experimental techniques and measurements have 
not reached a completely satisfactory level, describes a technical situation in need of an 
innovation. 

A new measuring technique based upon a cyclic perturbation of the incident radiation 
was initiated in the spring of 1966 at the NASA Lewis Research Center. A mathematical 
analysis of this new technique was first presented in reference 3. Previously obtained 
data using this method have been reported in references 4 to 8. The intent of this report 
is to summarize in a single, readily available publication the results obtained to date 
with the cyclic radiation technique. The present report therefore contains the following: 
(1) a complete analysis of the technique along with a discussion of its possible errors 
and its advantages and disadvantages, (2) an experimental validation of the method, and 
finally (3) experimental thermal radiative property data obtained with this method over 
a range of temperatures (150 to 600 K) and for a variety of materials. Metals, coatings, 
and a composite material (a CdS solar cell) were investigated, utilizing the cyclic tech- 
nique, to demonstrate its capabilities and its general utility. In addition, preliminary 
data are included on the specific heat of a paint as obtained by the cyclic radiation tech- 
nique. 
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THEORETICAL BASIS OF METHOD 


The cyclic method for measuring the normal absorptance and hemispherical emit- 
tance of surfaces is based on the temperature response of a thin, flat, isothermal model 
exposed to radiation (either solar or thermal) which is varied in a prescribed periodic 
manner after a steady-state isothermal model temperature T m has been established. 
The instantaneous heat-balance equation governing the model temperature is given by 

wc ^ - aA I(t) + q - 2eaA T 4 (1) 

p dt s s 


where q is the net heat absorbed (or lost) by the sample from all sources other than 
that from the incident radiation I(t) which can have a variety of functional forms such as 
step, ramp, or sinusoidal. (All symbols are defined in the appendix.) Now if it is spec- 
ified that the incident radiation vary in a convenient, periodic fashion with small ampli- 
tude about a mean level, the variation of the sample temperature about its equilibrium 
value T m will be small, and equation (1) can be linearized with respect to the equilib- 
rium temperature as 
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Since the model temperature is only perturbed, it is reasonable and expedient to assume, 
in each case, that the product of the temperature derivative and the delta temperature is 
small and can be neglected in the heat-balance equation. 

For a sinusoidal variation of the imposed radiant intensity 

I(t) = I Q + i sin u > t = I Q (1 + j sin wt) (3) 

equation (2) can be approximated as 

d(T - T ) 

wc p — 55- = oA s jI 0 sin u* - 8€crA s T^(T - T m ) (4) 
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Casting the original differential equation in this form illustrates one of the major assets 
of the cyclic approach for measuring thermal radiative properties. Most calorimetric 
methods must either assume that the heat-loss term in equation (1) is negligible or 
include a means to correct for it. In some instances the correction can be of the 
same magnitude as the terms of interest. In the cyclic technique the steady-state heat 
loss is unimportant because it is automatically accounted for by the mean temperature 
T m , which is an accurately measurable quantity. 

A sinusoidal variation of incident radiation is used since it is one of the simplest 
functions to handle analytically and, in addition, provides a convenient closed-form solu- 
tion. The solution of equation (4) then becomes 

T - T m = A sin(cut - (p) + Ce _t / 0 ' (5) 


where 
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Equation (5) is the complete solution of the differential equation. The constant C is de- 
termined from the initial condition: T = T m at t = 0. The last term of the solution de- 
cays exponentially with time and, after a sufficiently long interval, becomes negligible 
compared with the periodic term. Eventually, the sample temperature varies only in a 
steady-state sinusoidal fashion. Figure 1 illustrates graphically the variation of the 
model temperature on exposure to a sinusoidally varying incident radiation. It is noted 
that the sample temperature is out of phase with the incident radiation and lags behind it 
by an angle <p which varies from 0° at a low frequency to (but never exceeding) 90° for 
very high cyclic frequencies. The tangent of this angle is proportional to the cyclic fre- 
quency of the incident radiation as well as the time constant of the sample. 
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Figure 1. - Model temperature and incident radiation variations with time. 


The steady-state periodic solution is of interest because it offers a very effective 
and convenient means for determining simultaneously the radiation properties of mate- 
rials. Thus, by measuring the phase angle cp for a given frequency of incident radia- 
tion, the hemispherical emittance can be determined from 


coCpW 

8orA s T m tan V 


(6) 


and the absorptance may be determined from a measurement of the temperature ampli- 
tude and phase angle 


COCpW 

jl A 
J o s 


A 

sin cp 


( 7 ) 


Several observations concerning the analysis are appropriate. It should be recog- 
nized immediately that the accuracy of the thermal radiative property data is directly 
dependent upon accurate values for the specific heat of the material. Generally, accu- 
rate specific-heat data are available for pure metals; however, accurate values are not 
readily available for alloys, coatings, paints, and composites. 

The analysis assumes that the heat absorbed or lost by the sample q is independent 
of the sample temperature. This assumption must be looked at rather closely since it 
will affect the thermal radiative property determination especially at low temperature 
levels, where the amount of heat radiated away from the model is small. The concern 

here is the magnitude of the term (dq/dT) T (T - T ) in equation (2). K this term is 

m 

considered, the solution to the heat-balance equation is straightforward and similar to 
equation (5). The tangent of the phase angle is given by 
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tan cp = 
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where the plus or minus sign is for a heat gain or loss and K is the heat-exchange tem- 
perature coefficient K = (dq/dT) T 

m 

Solving for the true emittance yields 
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so that the difference between the measured emittance and the true emittance is given by 

e meas ~ e true ~ T 

8<A s T m 

It is noteworthy, however, that this difference is significantly less than that from the 
total heat-loss effect, which must be eliminated or accounted for in other methods for 
determining radiation properties. Furthermore, the effect of the heat-exchange temper- 
ature coefficient on the measured emittance can be minimized by using small-diameter, 
low-thermal-conductivity, low- emittance support wires or by using as large a sample 
surface area as possible. 

Unlike the emittance, the absorptance determination is not affected by the heat-loss 
temperature coefficient since the effects of a heat loss do not appear independently in 
equation (7) but are included in the temperature amplitude and the phase angle. As a re- 
sult the true absorptance is determined regardless of the heat-loss effect. 


UNCERTAINTY ANALYSIS 

An uncertainty analysis is useful to establish the capabilities and possible limita- 
tions of the cyclic technique and to indicate the major possible sources of error. In ad- 
dition, the instrumentation requirements for securing accurate data over a range of con- 
ditions are established. In the cyclic technique, as in all techniques, the desired pa- 
rameters (emittance and absorptance) are determined in terms of conveniently measur- 
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able parameters. The uncertainty analysis must therefore be made in terms of the esti- 
mated uncertainties inherent in these parameters. 

The uncertainty in the emittance and absorptance as determined from equations (6) 
and (7) is given by (ref. 9) 



Considering each term separately the uncertainty in w and A is due to the independ- 
ent measurements of the weight and surface area of the samples. For samples of rea- 
sonable size and with proper equipment, the sample weight and area can easily be deter- 
mined within an uncertainty of 0. 01. The uncertainty in the applied cyclic frequency 
5ct>/o> can also be kept to within 0. 01. 

The uncertainty in the specific heat of the material is difficult to estimate. For 
most pure metallic elements, data are available (ref. 10) which are accurate to within 
an uncertainty of 0. 02. However, for metal alloys, nonmetals, composites, or coatings, 
the uncertainty is considerably larger. For some materials, it may even be necessary 
to determine the specific heat, which for some materials may be obtained from the cy- 
clic technique (see later discussion). 

The uncertainty in the absolute temperature measurement is also difficult to esti- 
mate because of the thin samples (0.001 in.) and corresponding small- diameter thermo- 
couple wire (0.001 in.), both of which are often necessary. The combination of thin 
samples, small-diameter wires, and the many possible errors which can be introduced 
as a result of transient conditions precludes an accurate calibration of the thermocouple 
system. A best estimate of the uncertainty associated with the temperature measure- 
ment, based upon steady-state calibration of a typical system, is 2 K. 

The uncertainty in the temperature amplitude measurement 5A/A can be kept 
small because the measurement is relative and, with proper experimental technique, a 
high degree of resolution can be obtained over the range of the temperature perturbation. 
The uncertainty in temperature amplitude measurement is estimated to be within 0.01. 

Similarly, the uncertainty in the intensity amplitude factor can be kept quite small 
because the measurement is also relative and the resolution can be made quite large. 
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The uncertainty associated with this parameter is also estimated to be within 0.01. 

The final and most important parameter is the measurement of the phase angle be- 
tween the sample temperature and the imposed radiant intensity. The phase angle can 
be measured by a number of different methods; however, it is doubtful that this mea- 
surement can be made with an accuracy of better than 1°. The uncertainty associated 
with the determination of tan cp and sin cp required for an emittance and an absorp- 
tance determination is 


5 tan cp _ 25 cp 

tan cp sin 2 cp 

(13) 

5 sin cp _ 5 cp 

(14) 


sin cp tan cp 


Figure 2 shows the overall estimated uncertainty in emittance and absorptance as a 
function of the phase angle for temperatures of 100 and 300 K. The uncertainty in emit- 
tance is a minimum for phase angles between 30° and 60° but increases rapidly for 
smaller or larger phase angles and with decreasing sample temperature. For absorp- 
tance measurements, the experimental requirements are not as stringent. The uncer- 
tainty in absorptance is independent of temperature and decreases with increasing phase 
angle, becoming essentially independent of phase angle for angles greater than 45°. 
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Therefore, if only absorptance data are desired, the technique becomes relatively free 
of experimental difficulties and is ideally suited for absorptance measurements over a 
large temperature range. 


EXPERIMENT 

Experimental Requirements 

In order to effectively utilize the cyclic technique to provide thermal radiative prop- 
erty data over a range of temperatures, the experiment must be carefully established to 
satisfy the following primary requirements: 

(1) An environment which does not contaminate the sample surface and can be main- 

tained constant for long periods of time without introducing uncontrolled or var- 
iable heat- exchange effects between the sample and the environment 

(2) A sample mounting arrangement and sample temperature measurement which 

does not alter the isothermal sample assumption or nifluence the thermal re- 
sponse of the sample 

(3) A variable radiant intensity source which can be varied over a wide intensity 

range and also precisely controlled over a limited intensity range 

Figure 3 illustrates the requirements imposed by the sample response character- 
istics on an ideal experimental system for a range of temperatures. The two samples 
considered are 0. 001 -inch-thick aluminum and 0. 001-inch-thick aluminum with both 
surfaces painted with 3M Black Velvet paint. These two samples illustrate typical sam- 
ple response characteristics that can be encountered. 

For a practical intensity range from 0.01 to 10.0 solar constants (SC), the alumi- 
num sample temperature will vary from 145 to 690 K, while the painted sample will vary 
from 100 to 600 K. The corresponding time constants for aluminum (based upon Cp 
values from ref. 10 and emittance values from ref. 11) are 1000 to 8 seconds. For the 
painted sample (based upon Cp values from ref. 12 and an assumed emittance of unity), 
the time constants vary from 150 to 2. 3 seconds. The time constant is important be- 
cause it establishes the cyclic frequency required to maintain the proper phase angle and 
the minimum time required to complete an experiment. The cyclic frequency required 
to maintain a phase angle of 45° varies from 0. 001 radian per second (0. 57 rph) at 145 K 
to 0. 045 radian per second (26 rph) at 500 K for the aluminum sample, and from 0. 0143 
to 0. 27 radian per second (8 to 150 rph) for the painted sample at temperatures of 145 
and 500 K, respectively. The low cyclic frequencies and long time constants associated 
with the low sample temperatures require considerable flexibility and long-term stabil- 
ity in the experimental system. 
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Radiant intensity, 
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Figure 3. - Response characteristics of 0.001 -inch -thick aluminum and 0. 001 
inch-thick aluminum painted with 3M Black Velvet paint. 


Environmental Facility 

Figure 4 is a schematic drawing of the experimental facility being used. The test 
section is cylindrical with a 14- inch diameter and a 10-foot length. The lower 4-foot 
section of the facility is liquid- helium- cooled to provide cold wall pumping and a 4 K 
background environment to minimize the background radiant heat load to the samples. 
The vacuum capability of the liquid-nitrogen cryosorption roughing pumps and the liquid- 
helium cold wall is below 10 ^ torr. With such a vacuum, residual gas conduction is 
negligible. Furthermore, because no oil diffusion or roughing pumps are used, contam- 
ination of the sample surfaces by backstreaming of oil cannot occur. 

The sample test plane is located deep within the test section with liquid- helium- 
cooled baffles spaced along the length of the test section to minimize stray radiation 
from the uncooled upper sections of the facility. 
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Figure 4. - Schematic drawing of environmental facility. 
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Figure 5. - Schematic layout of sample mounting and thermocouple arrangement. 







Sample Mounting 


A typical mounting arrangement with four test samples is shown in figure 5. Gener- 
ally, four samples are installed in order that four sets of complete data may be obtained 
simultaneously. A radiometer is located in the same test plane at the center of the. 
sample array to monitor the radiant intensity level and to provide feedback control for 
the sinusoidal intensity variation. The samples are 1 centimeter by 2 centimeters and 
0.001 inch thick in order to achieve a fast response time. They are suspended by four 
0.001-inch Chromel or constantan wires spot welded to each corner of the samples. 

Two of the wires are a matched pair Chromel- constantan thermocouple. Chromel or 
constantan wire, both of which have low thermal conductivities, is used to minimize 
thermal conduction along the wire. A wire diameter of 0.001 inch is used to minimize 
thermal conduction and to minimize the effect of the wires on the thermal response of 
the samples. 


Solar Simulator 

The solar simulator consists of a 12-kilowatt carbon arc with associated optical 
elements (fig. 4). The simulator intensity output can be continuously varied by a mov- 
able zoom lens over the range from 100 to 250 milliwatts per square centimeter. Lower 
intensity levels can be obtained by using fine- wire- mesh screens that act as neutral den- 
sity filters (ref. 13). 

The sinusoidal intensity variation is provided by automatic control of the movable 
zoom lens. The control system is composed of a calibrated silicon solar cell which is 
coupled with a reference sine generator into a differential amplifier that actuates a bal- 
ancing motor to maintain the proper position of the zoom lens. The reference sine sig- 
nal can be varied to produce both the desired amplitude and the cyclic frequency of the 
sinusoidal perturbation. 

The spectral energy distribution of the solar simulator was measured with a double- 
pass spectrometer that had been calibrated by a transfer method over a wavelength 
range of 0.25 to 2.5 micrometers. That is, the spectrometer entrance slit was illumi- 
nated by a spectrally calibrated 1000-watt, quartz-iodine-tungsten-filament lamp whose 
calibration was traceable to a National Bureau of Standards (NBS) standard of spectral 
irradiance. The spectral distribution of the 12-kilowatt carbon arc solar simulator, ob- 
tained experimentally (ref. 14), is compared in figure 6 with the solar spectrum 
(ref. 15). Agreement of the carbon arc curve with the solar spectrum is good over the 
wavelength region of interest, 0. 35 to 2. 0 micrometers. 
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Figure 6. - Spectral irradiance of carbon arc solar simulator. 


Procedure 

The following procedure was used to obtain data: The radiant intensity was estab- 
lished at an appropriate level and maintained constant until the associated equilibrium 
sample temperature T m was reached. After equilibrium was achieved, the intensity 
was varied sinusoidally at a cyclic frequency and amplitude estimated as that required 
to provide the proper phase angle and temperature amplitude for accurate data reduction. 
(The four samples were selected to have compatible response characteristics so that 
data could be obtained simultaneously at a given cyclic frequency. ) The sinusoidal in- 
tensity perturbation was continued until cyclic equilibrium was reached as indicated by 
an unchanging temperature amplitude. Data for phase angle measurements were taken 
and the data reduced to determine if the phase angle was proper. If the phase angle was 
not suitable, the cyclic frequency and intensity amplitude were varied accordingly and 
data again taken at cyclic equilibrium. This procedure was repeated until the proper 
phase angle and temperature amplitude were obtained. 


Measurements 

The primary data recorded during the experiment was a Lissajous figure with the 
intensity on the x-axis and the sample temperature on the y-axis. A typical Lissajous 
figure is shown in figure 7. To obtain this figure a 0.001-inch-thick stainless- steel- 304 
sample was exposed to a radiant intensity having a sinusoidal radiant intensity perturba- 
tion factor of 0. 2 and cyclic frequency of 6 rph. The Lissajous figure provides all the 
necessary data on a single figure. The mean radiant intensity, intensity amplitude fac- 
tor, mean temperature, and temperature amplitude are all accurately recorded with 
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Radiometer output, mA 


Figure 7. -Typical Lissajous figure for 0. 001-inch-thick stainless- 
steel -304 sample. Cyclic frequency, 6 revolutions per hour; mean 
intensity, 2.59 milliamperes; intensity amplitude factor, 0.2. 


good resolution. The phase angle measurements are determined from the figure by 
(ref. 16) 


and 



( 15 ) 


sin cp = — 


( 16 ) 


The Lissajous figure also provides a useful monitoring function over the experiment 
in that the figure must close and repeat when cyclic equilibrium is achieved. 


Accuracy 


The accuracy of the emittance and absorptance determination can be estimated by 
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considering the uncertainty in the various required parameters that are determined from 
the Lissajous figures. The intensity and temperature level can be resolved quite accu- 
rately. Any errors involved in these parameters are due to the calibration of the radio- 
meter and the sample thermocouple. The uncertainty is within 0. 01 for the radiometer, 
and the uncertainty interval is within 2 K for the thermocouple. 

The uncertainties in tan cp and sin <p are determined from 


and 


(1 + tanV) 

tan (p 



(17) 


5 sin cp 
sin cp 



(18) 


The uncertainty in the intensity amplitude factors is estimated to be 6jI 0 /jI Q = 0.01 
and Si m /i m = 0. 02 for phase angles greater than 40°. For phase angles less than 40°, 
the uncertainty in 6i m /i m increases, while that for 6jl 0 /jl 0 remains constant. The 
uncertainty in the specific heat for stainless steel 304 is estimated to be 0. 05, while the 
uncertainty for a pure metal is estimated to be 0.02. Based upon these values, the un- 
certainties in emittance and absorptance as given by equations (11) and (12) are 
6e/e~0.09 and 5a/ a ~ 0.06 for the stainless-steel- 304 sample. However, had the 
sample been a pure metal, the uncertainties would have been 6e/e ~ 0.07 and 
5 a/ a ~ 0.04 because of the better estimate of the specific heat. 

In addition to the uncertainties which can be estimated, there are sources of error 
present that must be recognized and considered but whose effect cannot be estimated. 
Examples of these are nonuniformities in the radiant beam across the test plane and ra- 
diation bounce-back from the bottom of the tank to the samples. The uniformity of the 
beam has been evaluated in reference 14 and found to be within ±2 percent. The bottom 
and sides of the test facility have been painted black (a ~ 1.0), and apertures are used 
to minimize the bounce-back. 

Other possible sources of error are physical nonuniformity of the sample surface, 
contamination of the sample surface, and the spectral match of the solar simulator to 
the solar spectrum. These effects introduce systematic uncertainties into the radiative 
properties of the sample, and as a result, the radiative properties are only representa- 
tive for the source and the sample used. 
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RESULTS AND DISCUSSION 


Experimental Evaluation of Technique 

If the perturbation technique is to measure the total hemispherical emittance and 
normal solar absorptance accurately, the following requirements must be satisfied: 

(1) The phase angle, as determined from the Lissajous figure, should accurately 
represent the experimental system at thermal equilibrium. 

(2) Tan cp against u> (eq. (5)) and A against (jl 0 /w) sin cp (eq. (7)) should be linear 
relations for a constant sample temperature. 

(3) The heat-loss term q should be constant. 

(4) The associated radiation source should have good uniformity and very good sta- 
bility to ensure cyclic equilibrium at the lower sample temperatures. Since these prop- 
erties of the cyclic method are unknown, it was necessary to evaluate each of them 
under laboratory conditions to ascertain the validity of the cyclic radiation method for 
determining thermal radiative properties. 

In order to test the first property, a 1-kilowatt quartz- iodine-tungsten-filament 
lamp was used as a radiation source because of its very good stability over long periods 
of time. Its radiant intensity level was established to provide a sample temperature of 
approximately 260 K. After this temperature level was established, the radiant inten- 
sity was cycled at 4 rph with an intensity amplitude j of approximately 0. 565. The re- 
sulting Lissajous figure obtained for a 0. 001-inch-thick molybdenum sample is shown in 
figure 8. This figure was then analyzed and a phase angle cp determined. With this 
value for cp, a Lissajous figure was calculated and is shown as the dashed curve in fig- 
ure 8. Note that the agreement between the experimentally obtained Lissajous figure 
and the calculated figure is not exact over the complete cycle but the agreement is good 
and no asymmetries are evident in the experimental curve. The following observations 
can now be made: 

(1) Cyclic equilibrium was obtained since the figure closed and was repeatable. 

(2) The phase angle determined from one cycle accurately described the experi- 
mental conditions over one complete cycle. 

(3) The temperature response was truly sinusoidal. Thus, the linearized analytic 
solution is valid for the conditions of this particular experiment. 

(4) The resulting Lissajous figure has adequate resolution for an accurate determi- 
nation of the pertinent parameters T m , A, I i and jl 0 . 

Equations (6) and (7) reveal that tan cp and A must vary linearly with oo and 
(jI o /cu)sin cp, respectively, for the cyclic method to be valid. Experiments were con- 
ducted for cyclic frequencies varying from 2 to 30 rph, and intensity amplitude ratios 
ranging from 0. 0768 to 0. 565. These data are shown in figures 9 and 10 for 0. 001-inch 
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(a) Sample, 0.001-inch -thick molybdenum. Mean sam- (b) Sample, 0.005-inch -thick molybdenum painted with 

pie temperature, 275 K. 3m Black Velvet paint. Mean sample temperature, 

180 K. 

Figure 10. - Variation of temperature amplitude with cyclic frequency and intensity amplitude. Mean temperature 
held constant. 

molybdenum and for 0. 005-inch molybdenum painted with 3M Black Velvet. For every 
sample temperature, the data vary linearly for the materials used. Thus, this series of 
experiments demonstrates conclusively that the analysis, or the technique, is valid over 
a wide range of operating conditions. 

Although the analysis has been conclusively demonstrated to be proper, it remains 
to be proven that the emittance and absorptance determined are the true values. In par- 
ticular, the correctness of the emittance determination is directly related to the con- 
stancy of the heat-loss term q, as indicated in the theoretical section. If q varies 
slightly from a constant value, as would be expected for a temperature perturbation 
(e.g. , conduction losses), the true sample emittance as given previously by equation (10) 
is 


f — £ . 

meas true 


K 


( 10 ) 
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The effect of the variation of q with temperature on the measured emittance can be de- 
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Figure 11. - Variation of measured emittance with sample area. 


termined from a series of tests wherein the sample area A is the only parameter var- 

ied. With such data a plot of e „„ against the reciprocal of sample area is linear, 

moss 

and in particular, the intercept of such a plot is the true emittance, as shown in fig- 
ure 11. Such a series of tests were run at several temperatures, and the data obtained 
are presented in figure 12. The following observations can be made: 

(1) The measured emittance is related linearly to the reciprocal of the sample area 
as indicated by the analysis. 

(2) The slope of the straight line variation decreases as the mean temperature in- 
creases. Thus, the measured emittance and the true emittance are approximately equal 
at the higher temperatures. 

(3) At the lower temperatures, the true emittance and the measured emittance are 
different; therefore, a sample area as large as possible should be used. The maximum 
area usable is dictated, of course, by the uniformity of the radiation source (or the uni- 
formity of the sample temperature) since the sample must be isothermal for the analysis 
to be valid. 

It was stated previously in the theoretical section that the absorptance, as deter- 
mined by the cyclic technique, was unaffected by a heat-loss mechanism. In order to 
test this conclusion, the absorptance (or its equivalent, A/j sin <p) was also determined 
from the same series of tests used for figure 12(a). These data are plotted in fig- 
ure 12(b) with the reciprocal of the sample area as the independent variable. Each 
curve presented is for constant temperature, intensity and cyclic frequency. Therefore, 
the parameter A/j sin cp is directly proportional to the absorptance. For every tem- 


19 



II 


■ IIIII II i I ■ mil 


Mean cyclic 



°o 10 3 2 1 oo 10 3 2 1 

Sample area, A s , cm 2 


Figure 12. - Effect of sample area on measured parameters. 

perature, the absorptance is independent of the sample area and thus is independent of 
any heat-loss effect. 

Radiation source stability is a must since the material time constant increases rap- 
idly as the sample temperature decreases. In general, the time to achieve cyclic equi- 
librium is at least five times the material time constant. For example, a copper sample 
investigated at 162 K had a time constant of 1046 seconds. Thus, with a cyclic frequency 
of 0.0042 radian per second (2,4 rph), absolute control of the radiant intensity was re- 
quired for a period in excess of 2 hours to obtain a Lissajous figure that not only closed 
but also repeated itself. Time constants of this order do not represent a difficult exper- 
imental problem. However, for temperatures less than approximately 150 K where the 
time required to reach equilibrium is greater than 2 hours, minor intensity disturbances 
(or temperature disturbances) during a cyclic period can result in system perturbations 
which may require a restart of the test. 

Two methods have been used to minimize this problem. The time constant of a 
given material can be substantially reduced by depositing a thin layer of the material 
under investigation upon a metal substrate having a more suitable time constant. For 
example, molybdenum and lead can be used effectively because of their small time con- 
stants. If only emittance data are required, a tungsten- filament lamp, which is a very 
stable source, can be used at lower temperatures. 
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As a final check on the validity of the cyclic radiation technique, the solar absorp- 
tance of various materials was obtained from spectral reflectance data obtained with a 
double-pass monochromator and an 8- inch-diameter magnesium- oxide- coated integrat- 
ing sphere,, The integrating sphere method is a generally accepted method for deter- 
mining solar absorptance. The data obtained by using an integrating sphere are com- 
pared with those obtained with the cyclic technique in the following section. In general, 
the agreement between the data is very good for the materials chosen and the tempera- 
ture range covered. 

Although the data presented to establish the validity of the cyclic radiation technique 
for determining thermal radiative properties is not positively conclusive, the authors 
feel it is sufficient. As will be indicated herein, both accurate emittance and absorp- 
tance data can be obtained over a large temperature range for various materials, in- 
cluding metals, coatings, and composites, hi addition, the technique possesses a degree 
of versatility and convenience that no other method offers. 


Thermal Radiative Properties 

Metals . - The thermal radiative properties e, a, a/e of 12 metals (aluminum, 
copper, gold, molybdenum, nickel, platinum, silver, 304 stainless steel, tantalum, tin, 
titanium, and vanadium) were obtained. These data are shown in figure 13 for tempera- 
tures ranging from approximately 200 to 560 K. In general, the total hemispherical 
emittance of the metals increased with increasing temperature, the normal solar ab- 
sorptances were approximately constant, and the absorptance-emittance ratio decreased 
with increasing temperature because of the increase in emittance. 

In addition to the solar absorptance data obtained in the space environmental facility 
with the carbon arc solar simulator, solar absorptance data on the same metals were 
determined from spectral reflectance data obtained from integrating sphere measure- 
ments (ref. 8). The comparison between the measured carbon arc solar absorptance 
and the solar absorptance obtained from the integrating sphere reflectance measure- 
ments are shown in figure 13. For all the metals, the agreement between the absorp- 
tance values is good. 

Comparison of me tal data wi th theory . - There are three basic theories with which 
to conveniently compare experimentally determined total hemispherical emittance: (1) 
the Hagen-Rubens relation (HR), (2) the Drude single electron theory (DSE), and (3) the 
anomalous skin effect theory (ASE). A comparison of these theories with experiment is 
made in an effort to assess their suitability for making predictions of thermal radiative 
properties and to gain further insight into the validity of the cyclic technique. 

Each of the theories is based on the solutions of Maxwell's equation, which establish 
the correspondence between the optical parameters, refractive index n and extinction 
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coefficient k, and the electric and magnetic properties of a material. From the Fresnel 
equation for reflectance, the spectral normal emittance in terms of the optical param- 
eters is 


e (X) = (19) 

n 9 2 

(n + ir + k 

The total normal emittance is then determined by the integration of the product of 
the spectral emittance and the Planck spectral energy distribution over all wavelengths. 
Finally, the total hemispherical emittance for metals is determined from the total nor- 
mal emittance by (ref. 17) 


e = 1. 33 e n (20) 

A detailed discussion of the basic theories and their application to total emittance 
calculations is presented in references 18 and 19. 

The total hemispherical emittance data obtained experimentally for all the metals 
are compared in figure 14 with the analytical relation derived from the Hagen-Rubens 
spectral emittance. The Hagen-Rubens equation for the total hemispherical emittance 
is (ref. 18) 


e 



(21) 



Figure 14. - Comparison of ail experimental data with 
Hagen-Rubens equation. 
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The Hagen-Rubens equation correlates the absolute level of the total hemispherical 
emittance data for all the metals except titanium to within ±50 percent. However, it 
should be recognized that titanium is a very active metal so that surface contamination 
is probably the cause of the large disagreement noted for this metal. (The values used 
for direct- current electrical conductivity in equation (21) were obtained from refs. 10 
and 11 except for tin and 304 stainless steel, which were obtained from ref. 20.) For 
many engineering purposes, the comparison shown in figure 14 is both quantitatively and 
qualitatively acceptable. The fact that the theory does correlate the data for the 12 dif- 
ferent metals with modest success probably indicates that the direct- current electrical 
conductivity of the metal represents the effective conductivity of most metals regardless 
of the atomic structure and the small amount of surface contamination present. The 
comparison, however, leaves much room for improvement if very accurate thermal 
emittance data are required. In that case, experimental data on the particular material 
are needed. 

The experimental data for four metals are compared with the Hagen-Rubens, DSE, 
and ASE theories in figure 15. The metals chosen are the noble metals copper, gold, 
and silver and the nonnoble metal aluminum. The noble metals were chosen because of 
their compatability with the assumption of a single free conduction electron in the DSE 
and ASE theories (ref. 21). Aluminum was chosen because of its general application in 
many areas of space technology. 

Several observations can be made from the comparisons in figure 15. The Hagen- 

Theory 

• Hagen-Rubens (HR) 

Anomalous skin effect (ASE) 

Drude single electron (DSE) 



Sample temperature. K 
(c) Silver. (d) Aluminum. 

Figure 15. - Comparison of experimental data with theory. 
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Rubens theory predicts the absolute level of the emittance for copper and aluminum with 
fair accuracy but only approximately for gold and silver. The Hagen- Rubens prediction 
agrees more closely with the experimental data as the temperature approaches 250 K. 
Both the ASE and DSE theories predict emittance values which are lower than the exper- 
imental data. The variation of emittance with temperature for aluminum is predicted 
well by the Hagen-Rubens relation, while the ASE and DSE theories predict the slope 
much better for the noble metals.. 

Although the DSE and ASE theories are considered more rigorous than the HR 
theory, the quantitative agreement with experimental data is not as good. However, this 
disagreement must be viewed in light of the assumptions made in deriving the electro- 
magnetic theories. All the theories assume an ideal interaction between the incident 
electromagnetic waves and the material surface. Physically this means, for example, 
that the surface is optically smooth and uncontaminated. The fact that ideal surfaces do 
not exist in experiments is probably responsible for the differences noted between theory 
and experimental data. Thus, such factors as impurities, contamination, surface finish, 
and crystal structure can significantly affect the actual thermal radiative properties. 

Such nonideal surface effects tend to increase the emittance of metals. For tempera- 
tures greater than approximately 200 K, the HR approximation yields results higher than 
those of the more exact theories. Therefore, the HR relation frequently correlates the 
experimental data better than the more exact theories. This improved agreement is 
therefore probably fortuitous, but it is representative of results obtained with real ma- 
terials. 

In summary, all the metal data obtained by the cyclic technique have indicated the 
following conclusions: When compared with theory, the level of the data is proper, and 
in particular, the trend of the data with temperature is proper. In addition, the absorp- 
tance data obtained for 12 metals by utilizing the cyclic technique is in good agreement 
with data obtained from spectral reflectance measurements made using an integrating 
sphere. These results indicate that the cyclic method of determining thermal radiative 
properties is a valid one. 

Comp osite - CdS solar cell . - The optical properties of a composite material, 
namely a CdS solar cell, were obtained at temperatures corresponding to radiant inten- 
sities that varied from 0. 04 to 1. 1 SC. The time constants (eq. (5c) as determined from 
the measured phase angle cp and imposed frequency co are shown in figure 16. It is to 
be noted that the response time increases with decreasing temperature, and that the 
longest time constant obtained is 143 seconds at a cell temperature of 156 K. Therefore, 
the longest time required to reach cyclic equilibrium is of the order of 12 minutes. 

Since this time is relatively short, there are no difficulties in meeting the critical re- 
quirement that all experimental parameters be held constant for this period. 

Two data points obtained with the 1-kilowatt quarts- iodine-tungsten-filament bulb 
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Cell temperature, K Cell temperature, K 

Figure 16. - Time constant of cadmium sulfide solar cell. Figure 17. - Emittance of cadmium sulfide solar cell. 

O- 



Cell temperature, K Cell temperature, K 

Figure 18. - Solar absorptance of cadmium sulfide solar cell Figure 19. - Absorptance-emittance ratio for cadmium sulfide 

(carbon arc source). solar cell. 


are also included in figure 16 since the cell time constant is unaffected by the source 
spectral distribution. Data from the 1-kilowatt bulb are much more accurate than those 
obtained with the arc lamp because there are no instability problems. Despite this, 
however, all the time constants are very consistent regardless of the source used. 

With the time constant established, the emittance can be determined (eq. (6)) from 
the same data points. Again, data from both sources can be used to obtain emittance. / 

Solar absorptances are obtained by substituting in equation (7) the additional measure- 
ments of temperature amplitude A and intensity amplitude jI Q . In this case, however, 
the 1 -kilowatt- bulb data cannot be used since its spectral distribution is not comparable 
to the solar spectrum. 

The emittance and absorptance data obtained are presented in figures 17 and 18 in 
terms of the cell heat capacity per unit area. Unfortunately, this had to be done because 
no specific-heat data are available for CdS solar cells. As a result, the character of 
these plots is influenced greatly by the cell specific heat variation with temperature. 
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Even so, the plots could be utilized to obtain absorptance and emittance by those who are 
willing to estimate the specific heat of the solar-cell materials. 

The important optical parameter, the parameter establishing the cells operating 
temperature, is the absorptance-emittance ratio. This parameter can readily be ob- 
tained, without any assumptions, by ratioing the curves presented in figures 17 and 18. 
The absorptance-emittance ratio is shown in figure 19 for cell temperatures ranging 
from approximately 150 to 325 K. The important trend to note is the rapid increase in 
the absorptance-emittance ratio as the temperature is decreased. As the temperature 
increases from 155 to about 200 K, the absorptance-emittance ratio decreases from 0. 70 
to about 0.43, a decrease of 39 percent. A further increase in temperature produces 
little change; the a/e ratio is approximately constant at a value of 0. 43 from 200 to 
325 K. Comparison of the 300 K experimental point presented in reference 22 with the 
preceding data shows good agreement. As a point of information, it should be noted that 
the cell was not producing power when the optical properties were obtained. As a result, 
the temperatures measured are slightly high. Calculations indicate that, for a cell hav- 
ing an electrical conversion efficiency of 5 percent, the equilibrium cell temperature, 
while producing power, may be as much as 6 K lower than that presented herein. 

Paints and coatings . - Currently underway at Lewis is a program to determine the 
thermal radiative properties of various thermal control paints and coatings using the 
cyclic radiation technique. The cyclic technique can be used readily if the coating (or 
paint) can be properly applied to a suitable substrate and if the specific heat of the paint 
or coating is known. This latter requirement presents a difficulty since the specific 
heats of paints and coatings are not generally available. Fortunately, however, the cy- 
clic technique can be used in several ways to determine the specific heats of paints and 
coatings, so that thermal radiative properties can be determined. To illustrate the use 
of the cyclic technique for coatings and paints, data were obtained for 3M Black Velvet 
paint and for the Z-93 white coating. The specific heat for Z-93, which was unknown, 
was also determined in this experiment. These data, along with the solar absorptance 
data obtained at 300 K by integrating sphere reflectance measurements, are presented 
in figure 20. 

3M Black Velvet paint: The 3M Black Velvet paint sample consisted of a 0.005-inch- 
thick molybdenum substrate upon which the paint was applied to both surfaces. The 
paint weight was approximately 8 percent of the substrate weight. The specific heat of 
the paint was obtained from reference 12. The radiative property data indicate that the 
emittance and absorptance are approximately equal and are greater than 0. 95 over the 
temperature range from 180 to 315 K. The comparison of the absorptance determined 
by the cyclic method with the integrating sphere reflectance measurement (a = 0. 97) at 
a temperature of 300 K is good. The measured cyclic absorptance of 1.06 at 315 K is, 
of course, too high. This value, however, is within the accuracy of the method and in- 
dicates that the solar absorptance of the 3M Black Velvet paint is very nearly unity. 
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Figure 20. - Absorptance and emittance of 3M Black Velvet paint 
Z-93 white coating. 


Z-93 white coating: The method chosen initially to determine the Z-93 white coating 
specific heat was as follows. The coating was applied to both surfaces of two different 
substrate materials. The substrate materials selected were 0. 005- inch-thick nickel and 
molybdenum and were chosen because of their different thermal response characteristics. 
The specific heat of the coating was determined by a simultaneous solution of the meas- 
ured sample time constants for a given temperature and for an assumed equal coating 
emittance on the samples by 



( 22 ) 


The data in figure 20 show the emittance of the Z-93 coating to be approximately 
0.92 and the absorptance to be 0.23 over the temperature range from 165 to 230 K. The 
comparison between the absorptance data from the cyclic technique (which is the average 
of the two samples) and the integrating sphere reflectance measurements at a tempera- 
ture of 300 K (a - 0. 18) was not as good as expected. The differences noted here are 
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attributed to differences between the arc lamp spectral distribution and the solar spec- 
trum (ref. 14), errors in the wavelength limitations of the integrating sphere instrumen- 
tation (see ref. 8), and the error in the specific-heat determination. 

The determination of the coating specific heat is another illustration of the versatil- 
ity and usefulness of the cyclic technique. However, it should be noted that the proce- 
dure of using two different substrates to determine specific heat is not the only proce- 
dure that can be used. In fact, it probably is not the best procedure because the as- 
sumption of equal coating emittance of the two samples is subject to question. A more 
suitable approach is to coat the back side of a pure metal standard whose absorptance is 
known and expose it to a known cyclic radiation. The coating specific heat can then be 
determined from 



a s A s n q sm^ 

w c u>A 



(23) 


where a„ is the solar absorptance of the metal standard. Thus, the coating specific 
heat can be determined in terms of the absorptance of a pure metal, which for most 
metals is independent of temperature. 


CONCLUDING REMARKS 

Considerable experience has been gained in the application of the cyclic radiation 
technique to determine thermal radiative property data for various materials. The pri- 
mary advantages of the method are 

1. The emittance and absorptance are determined simultaneously. 

2. Radiative property data can be obtained conveniently for many different types of 
materials over the temperature range from 150 to 600 K. 

3. The technique is not affected by the magnitude of any extraneous heat-loss ef- 
fects. The variation of the heat loss over the range of the temperature perturbation 
does not affect the absorptance determination but can introduce a slight effect into the 
emittance determination. The error in emittance, however, is much less than that en- 
countered by other methods, which must account for the magnitude of the heat loss. 

4. There is versatility in the technique since the experimental parameters, cyclic 
frequency and intensity level, can be controlled to yield accurate data irrespective of the 
temperature level and the magnitude of the material radiation properties. 

5. The specific heat of coatings and paints may be determined by using a suitable 
metal standard for a substrate. 
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The cyclic method is fairly convenient to apply but there are certain experimental 
requirements that must be met, as follows: 

1 0 The sample configuration is governed by the requirements of the technique. The 
sample must be a thin, flat section for reasonable response time, and the sample must 
have an exposed area of at least 2 square centimeters. 

2. The sample material (or the substrate) must be suitable for attaching small- 
diameter support and thermocouple wires. 

3. The specific heat of the sample material (and substrate) must be accurately 
known. 

4. A versatile radiation source must be available which can be precisely controlled 
in a sinusoidal manner. In addition, the amplitude and cyclic frequency of the sinusoidal 
variation must be variable over a large range and capable of being maintained constant 
for long periods of time. 

Considerable data have been obtained to check the validity of the method and no 
serious problem areas have arisen. The method has generally provided very consistent 
data over a large range of operating conditions and experimental variables. The pri- 
mary requirement for obtaining accurate data is to maintain the phase angle between 40° 
and 60°. This is conveniently accomplished by adjusting the incident radiation cyclic 
frequency. 

The most serious limitation of the method is considered to be the fairly long re- 
sponse times encountered for some materials, during which absolute control over the 
experiment must be maintained. Therefore, thin samples having the desired fast re- 
sponse times are preferred and, in general, have been used. 

Thermal radiative property data were obtained for a variety of materials at temper- 
atures ranging from 150 to 600 K. The following summary remarks can be made. 

The total hemispherical emittances and solar absorptances obtained for 12 metals 
are consistent and follow the accepted variation with temperature. No comparison can 
be made regarding the level of the emittance obtained except that the results appear to 
be reasonable when compared to theory and other available experimental data. The com- 
parison between the cyclic solar absorptances and those obtained from integrating sphere 
reflectance measurements is good. 

A composite, namely a thin film CdS solar cell, was investigated over temperatures 
ranging from 155 to 325 K. The absorptance-emittance ratio a/(e^ + e b ) of the cell was 
approximately constant over the temperature range 200 to 325 K. Below 200 K the ratio 
increased rapidly. The rapid increase was attributed to a decrease in emittance below 
a sample temperature of 200 K. 
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Preliminary radiative property data were obtained for a Z-93 white coating and a 
3M Black Velvet paint. In addition the cyclic technique was used to determine the spe- 
cific heat of the Z-93 coating. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 17, 1969, 

124-09. 
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APPENDIX 


SYMBOLS 


A amplitude of temperature perturba- 
tion during cyclic period, K 

A g sample area, cm 2 

C constant of integration 

c specific heat, J/(g)(K) 

Ir 

I(t) time varying radiant intensity, 
J/(cm 2 )(sec) 

2 

I mean radiant intensity, J/(cm )(sec) 

i amplitude of intensity perturbation, 
J/(cm 2 )(sec) 

i intensity amplitude at mean temper- 
ature during cyclic equilibrium, 
J/(cm 2 )(sec) 

j intensity perturbation factor, i/I 

K heat-exchange temperature coef- 

ficient, J/(sec)(K) 

k extinction coefficient 

n refractive index 

q all heat- exchange terms other than 

those specified by remaining 
terms in eq. (1), J/sec 

T test sample temperature, K 

T time rate of temperature change, 

K/sec 

T m mean cyclic temperature of sample, 

K 


t time, sec 

w test sample weight, g 

a solar absorptance 

6( ) uncertainty interval 

e total hemispherical emittance of 

sample 

e n total normal emittance 

6' time constant of material, sec 

X wavelength 

a Stefan- Boltzmann constant, 

J/(cm 2 )(sec)(K^) 

a electrical conductivity, (ohm- cm) - * 

cp phase angle between sample tem- 

perature and cyclic incident 
radiation, rad 

a> cyclic frequency, rad/sec or revo- 

lutions per hour (rph) 

Subscripts: 

b back surface 

c coating 

f front surface 

m mean temperature 

s standard 

1 substrate 1 

2 substrate 2 
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